Nonregenerative MIMO Relaying with Optimal 
Transmit Antenna Selection 



Steven W. Peters, Student Member, IEEE, and Robert W. Heath, Jr., Senior Member, IEEE 



Abstract 

We derive optimal SNR-based transmit antenna selection rules at the source and relay for the 
nonregenerative half duplex MIMO relay channel. While antenna selection is a suboptimal form of 
beamforming, it has the advantage that the optimization is tractable and can be implemented with only 
a few bits of feedback from the destination to the source and relay. We compare the bit error rate of 
optimal antenna selection at both the source and relay to other proposed beamforming techniques and 
propose methods for performing the necessary limited feedback. 

I. Introduction 

Despite the lack of precise knowledge of its basic theoretical behavior and limits, relaying is beginning 
to find practical application in standards such as IEEE 802. 16j [1]. By deploying relatively inexpensive 
relays, service providers can reduce the number of base stations required to serve a given area, or increase 
capacity at the cell edge. 

Relaying research efforts have also increased recently [2]-[7]. Capacity bounds for the full-duplex 
MIMO relay channel were derived in [2], [3]. The authors of [6] derive the optimal infinite-SNR 
diversity-multiplexing tradeoff for the half duplex MIMO relay channel and find that a compress-and- 
forward strategy is optimal in this sense. Recently, practical strategies have been developed for MIMO 
relaying. Both [4] and [7] derive the mutual-information-maximizing nonregenerative linear relay for 
spatial multiplexing when the direct link is ignored. 

This letter derives the optimal transmit antenna selection criteria at both source and relay; i.e., all 
transmissions occur using the transmit antenna that will give the destination the highest post-processing 
signal-to-noise ratio. We consider the case where only a single spatial stream is to be sent from source 
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to destination. This scenario arises when the channel is ill-conditioned (i.e., there is a dominant path of 
propagation in the source-destination channel), or if robustness via diversity is preferred over throughput 
(i.e., near the cell edge). 

UnUke most previous practial MIMO relay results (e.g., [4], [5], [7]), the strategy derived here is the 
optimal transmit anteima selection strategy when the direct Unk from source to relay is not ignored. We 
prove that transmit antenna selection, combined with an MMSE receiver at the destination, achieves the 
full diversity order of the MIMO single relay channel. That is, at high SNR the probability of outage 
decays with SNR as quickly as is possible in such a model. Further, antenna selection requires less 
feedback than beamforming. Distributed space-time codes, which may also achieve the full diversity 
gain, not only require their own level of overhead for coordination and synchronization, but also require 
the relay to be able to decode the message transmitted by the source. 

Compared to recent results using limited feedback beamforming [8], under the tested parameters given 
in the aforementioned paper, antenna selection at both source and destination is about twice as likely to 
cause bit errors as a Grrassmaimian codebook with 16 codes, which is a loss of about 1 dB at high SNR. 
In return, anteima selection requires only log2 NsNr bits of feedback versus 3 log2 N + 2b bits in [8], 
where Ns and Nr are the number of antennas at the source and relay, respectively, N is the size of the 
Grrassmaimian codebook, and b is the quantization in bits of the SNR feedback required in [8]. 

This letter uses capital boldface letters to refer to matrices and lowercase boldface letters for column 
vectors. The notation ||h|| refers to the L2-norm of the vector h, and H* is the complex conjugate 
transpose of the matrix H. The vector h^*) refers to the ith column of the matrix H. Finally, 

A = B^ lim i^ = -5. 

SNR^oo log SNR 

II. System Model & Antenna Selection 

We assume a single source S transmitting information to a destination V with a single relay TZ 
aiding the transmission. The source, destination, and relay are equipped with Ns, Nd, and antennas, 
respectively. All nodes operate in half-duplex mode. Unlike most prior work in MIMO relaying, we do 
not ignore the direct link between S and V. 

The source S wishes to transmit the scalar symbol s to V, where E|sp = Eg = SNR, Eg is the 
average power constraint at both S and TZ, and cr^ = 1 is the overall noise power at each node. Since 
the signal-to-noise ratio is the metric of interest, an imbalance of noise energy among the nodes can 
be modeled in the appropriate fading parameter for Hxy- For instance, if the relay has noise power 
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(T^, in an independent Rayleigh fading environment these definitions would change the channel fading 
parameter of the corresponding exponential distribution from Xsr to Xsrcj'^. 

We denote the channel from X Xoy, X e {8,11}, y G {7^,P}, A" / 3^, as Hxy, and h^y is the 
vector channel from the zth transmit anteima at X to y. We also define 



(.) - 2 



TXY 



SNR (1) 



to be the equivalent receive SNR from Xi — > y. 

We assume the block fading model. In the first stage, if S transmits s on anteima i, TZ receives the 
signal 

where ur is the zero-mean spatially white complex Gaussian noise vector with covariance ct^Inr as 
observed by TZ. Since the relay is also transmitting on only one of its antennas, it must combine its 
received vector to form a single symbol. It can be shown that the optimal way to do this is to perform 
MRC on the signal, resulting in a scalar 

sr = Q!(hl!]^)*yR, (3) 
where a is the scaUng factor to ensure TZ transmits at its expected power constraint; i.e., 

= —T\ ^-T^ • (4) 

l|hgill^ + l|h«||VSNR 

At V, the first stage results in 

yD,i = h%s + nD,i. (5) 
In the second stage, TZ transmits sr to V on anteima k: 

YD, 2 = ^rdSR + n£>,2. (6) 

The destination now has two observations containing s. To put the channel in standard MIMO notation, 
we define 

h = I 1 (7) 

Vllh^-iP+l/SNR 



l|hL'illVll4'ilP+l/SNR 

yn = I '''' ] (9) 

V yD,2 J 
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SO that 

yD = hs + n. (10) 

We assume the destination V now applies a linear filter w to y/j to obtain an estimate of s. Although 
suboptimal, we will see later that in some cases the destination may wish to apply MRC (w = h) on 
Yd, and doing so would result in the post-processing signal-to-noise ratio 7^*^ of ([TT]) at the top of the 
page. In this form, it is easy to see that, if 

{tE < ill} n {^'rI > 7fi(7a + mig - ig)} , (12) 

then 7^*) < 7^*]-, and relaying is worsening performance. This occurs when the SNR from 7^ to 2? is very 
good relative to the others, and the SNR from 5 to 7^ is worse than the direct SNR. Effectively, the TZ to 
V channel is dominating the received signal, but it consists of mostly noise relative to the direct signal. 
Recall that MRC is only optimal when the observations contain the same noise variance [9]. Because of 
the amplified noise at TZ, this is not the case here. In this case, one can show that the optimal receive 
filter in the minimum mean-squared error (MMSE) sense is 

w = R-^Ry,„ (13) 

where Ry„ = Elyoy})} and Ry^s = Ejy/js*}. The post-processing SNR is then 

(i) (k) 

7 -1SD+ (i) (k) ^^^^ 
TSR + 7Rf) + 1 

Note that this requires the destination to have knowledge of ||hs/j||. If this is not possible, suboptimal 
MRC resulting in the SNR of ([iTI ) may be used instead, which requires less training. A method for 
obtaining this CSI is presented in Section |llll 

Note that in ([141 ). for fixed jgj^ and Jgl^, 7^'^ is maximized when 7^^^ is maximized. Thus, the antenna 
selection at the relay is independent of the selection at the source, and we can substitute the index of the 
optimal relay transmit antenna ko in for k in all subsequent equations. The same cannot be said of the 
regular MRC equation (fTTl ). 

Finally, we note that antenna selection at the relay is suboptimal, and the optimal strategy in this case 
is intuitive; since the SNR expression (fT4l ) is the addition of the independent SNR terms for the parallel 
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channels to the destination from the source, the relay should apply a filter that maximizes the SNR to the 
destination. One can show that this filter is W = v(^)(h^^^)*, where v^^^ is the right singular vector of 
Urd corresponding to its largest singular value, and io is the index of the source antenna that maximizes 
([T4l ). Intuitively, W is the combination of a receive filter matched to Hsr and a transmit beamforming 
vector matched to Urd- Implementing this filter would require perfect knowledge of Urd at the relay 
and an SVD operation. All of our results hold with this optimal strategy, with ||h^j^^|p replaced with 
^RD = o'rd' square of the largest singular value of Urd- 

III. Training and Limited Feedback 

We now discuss how channel state information might be obtained in the channel of interest so that 
a reliable antenna selection strategy may be implemented. All three channels need to be estimated at 
their respective receivers; this can be accomplished using previously studied MIMO training methods. 

(i) 

Only knowledge of the link SNRs (i.e., 7)^y's) is required for transmit antenna selection. Therefore a 
low complexity signal, such as a short narrowband tone, may be used for estimating SNR to choose 
an antenna to train from. This is first sent from 7^ to P from each relay antenna. V then feeds back 
which antenna TZ should use to transmit, and, from this antenna, a training sequence suitable for channel 
estimation is sent to the destination. The source repeats this process with its transmit antennas, with the 
relay forwarding its received signal on its optimal antenna. This way, the destination can estimate the 
SNR between the source and relay to perform MMSE combination as described earlier. 

The destination finds (fT4b for each source antenna, then feeds back to the source the index of the 
antenna that resulted in the largest 7^*). The source then transmits a training sequence from this antenna, 
which does not need to be forwarded by the relay. This process requires log{NsN£,) bits of feedback, 
two time slots of training, and Nr + 2Ns time slots for SNR estimation. Minimizing the time required 
for SNR estimation is thus important for this feedback strategy. 

IV. Diversity Analysis 

Antenna selection is used to exploit the diversity gain available in the channel. Using ([141 ) we now 
show that this strategy achieves full diversity gain. We first give an upper bound on the diversity order of 
the half-duplex MIMO relay channel when the source and destination transmit orthogonally in equal time 
slots. Yuksel and Erkip [6] have derived this result for arbitrary time sharing when the source is allowed 
to transmit in the second time slot, so this result is a special case of their derivation. This derivation 
is included here to prove that our added restrictions (i.e., equal transmission times, source silent in the 
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second time slot) do not decrease the maximum diversity order of the channel. We first define 



Ibc = I{S;Yr,Yd,i) (15) 

h = I{S;Yd,i) (16) 

h = I{Sr;Yd,2) (17) 

Imac = h+h, (18) 



where S is the random variable corresponding to the transmitted signal from the source, Yr is the 
received signal at the relay, Yo.n is the received signal at the destination in the nth time slot, and Sr is 
the transmitted signal at the relay. Using equations (27) and (28) in [6] with t = 0.5 and the source not 
transmitting in the second time slot, 

1(3; Yd) < 0.5mm{lBc, hiAc}- (19) 
Now we can bound the probability of outage for a fixed Iq as 

Pout = MHS;Yd) < lo} 

> Pr {0.5 mm{lBC, Imac} < h} ■ (20) 

The event where the minimum of two variables is less than a constant is equivalent to the union of 
the events that each of the variables is less than the constant. Defining Pout,BC = P^{Ibc < 2/o)» and 
similarly for Pout,MAC^ we can write 

Pout > Pr({/Bc < 2/o}|JUmac < 2/o}) (21) 

= Pout,BC + Pout,MAC — 

Pr (^{Ibc < 2/0} f]{lMAC < 2/0}) . (22) 

Recall from (fTSl ) that Imac is the sum of two nonnegative random variables. Such a sum is always less 
than or equal to twice the maximum of the two random variables. Then, by making the codebook for 
Sr independent from that of S, and defining Poutji — < Iq) and Pout.i^ similarly. 

Pout > Pout,BC + Pr{max{/i, 12} < Iq} - 

Pr ({/sc < 2/0} f]{lMAC < 2/0}) (23) 

> Pout,BC + Pout,hPout,l2 ~ 

Pr ({Ibc < 2/0} f]{lMAc < 2/0}) . (24) 
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Conversely, the sum of Ii and I2 is always greater than the maximum of the two. Also, note from (ITSl) 
and (O that Ibc > h so that 



Pout > Pout,BC + Pout,IiPout,h ~ 

Pr [{Ibc < 2/o}f|{max{/i,/2} < 2/o}^ 

= Pout,BC + Pout,hPout,l2 ~ 

Pr{max{/i3C,/2} <2/o}. (25) 
Finally, again assuming independent channels on all links. 

Pout > Pout,BC + Pout,hPout,h ~ 

Pout,BC^r{h < 2/0}. (26) 
From MIMO information theory we know that (see [6], Sec. Ill and IV) 

Pv{Ibc<c} = Ns{Nr + Nd) (27) 
Pr{/i < c} = NsNd (28) 
Pr{/2<c} = NrNd, (29) 



for all c G M. Thus, the last term in (126) will decay as NsNr + NsNd + NrNd with log SNR and is 
thus irrelevant to the diversity analysis. The first term will decay as Ns{Nr + Nd), while the second 
term decays as N£){Ns + Nr), so that 

Pout < NsND + NRmm{Ns,ND}. (30) 

We now derive a lower bound on the diversity order of optimal antenna selection in flat i.i.d. Rayleigh 
fading by using (fT4l) . First we define 

lM^^=llWRh^/illUlt^ + ^)- (31) 
Since we choose the source transmit antenna that maximizes the SNR 7 at the destination. 

Pout = Pr{7 < 70} 

= Pr{max{7*^*)} < 70} 

i 

= Pr{max{7ji+7j],'^)}<7o}. (32) 
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As before, the sum of two random variables is greater than the maximum of the two. 

Pout < Pr{max{max{7^*]5,7^*]^£,^}} < 70} 

= FT{m^x{jg,^l'^^}<^o}- (33) 

Since each channel is mutually independent of the others, and the channel from each source antenna to 
the destination is also independent from the others, we define Pout.SD = ^"^{l^sD ^ ^^)' ^^^^ 

Pout < Pr{max{7j]5} < 7o}Pr{max{7^*]j5^} < 70} 

= (Po«t,5D)'^"Pr{max{7j]j'^£,)} <7o}. (34) 

Now define 

7M,i = min{7f jj , 7 J^,^ } . (35) 
If 7A/,i > 1, then ■^^g^^ > 7M,i/3. Otherwise, Jg'^^ > (7M,j)^/3. In either case, since 70 is arbitrary. 



i 



we let 70 > 1/3 and proceec 

Pout < (Po«t,5D)^" Pr{max{7M,J < 370}. (36) 
We can again split up the minimum event into a union: 



Pout < {Pout,SDf' X 



Pr({max{7« } < 3^o}[j{rnf^{l^Rh^} < 37o}) 

{Pout,SDf' X 

Pr{max{7« } < 370} + Pr{7&^ < 37o} - 
Pr{max{7« } < 37o}Pr{7iy < 37o} • 



(37) 



Again, since the channels between each source transmit antenna and the relay are independent, we define 

Pout,SR = Pr(7sR < 370) and Pout,RD = ^^{Ird < 37o), and 

Pout < {Pout,SDf' X 

iPout,SR)^^ + Pout,RD — 



{Pout,SR)'^^ Pout,RD 



(38) 



'since Pout is monotone increasing witli increasing 70, no loss in generality occurs by assuming 70 > 1/3. For example, let 
7L < 1/3. Then Pr(7 < jl) < Pr(7 < 70). Thus, if Pout = d, then Pr(7 < 7L)>d. 
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where again the last term will decay much quicker than the others and can be ignored. The first term, 
after multiplication, will decay as NsNo + NsNr, while the second term decays as NsNd + Nj^Nd. 
Thus, 

Pout > NsND + NRmm{Ns,ND}. (39) 

Combining ( [39l ) and (l30l ) we see that the proposed antenna selection achieves the full diversity gain in 
the channel. 

V. Simulation Results 

We present a simple simulation to compare to a recent result on limited feedback beamforming [8]. 
For each case shown, we simulate the relay channel with Ns = = Nd = 3 using BPSK modulation 
and an i.i.d. Rayleigh channel at each link. Bit error rate (BER) is the metric of interest. Figure |2] gives 
the results for E{7j]j} = ^{Ird} = 2 dB for various £{7^*15}. Note that this graph corresponds exactly 
to Fig. 9 in [8], and we have included their results for a Grassmannian codebook with more than 20 bits 
of feedback. Using antenna selection at both S and TZ requires 4 bits in this case and results in a loss 
of approximately 1 dB at high SNR. 

The theoretical lower bound of Figure |2] is when the source can simultaneously beamform the BPSK 
symbols to both the relay and destination; obviously this is an impossible task. The "optimal" performance 
curve was found numerically in [8] using gradient descent to find a local optimum. 

Figure [3] shows the BER of uncoded BPSK versus Es/Nq for a relay channel with two antennas at 
each node. Note that increasing Es/Nq implies an increase in SNR at each link (recall that noise terms 
are normalized and E|sp = Es = SNR). The figure was generated using Monte Carlo simulations using 
10^ channel realizations for accuracy at high SNR, and demonstrates that antenna selection achieves the 
maximum diversity order available in the channel. 

VI. Conclusion 

We explored antenna selection as a practical way of achieving the full diversity order of the nonregen- 
erative MIMO relay channel. It was shown to achieve this diversity with a small SNR penalty relative 
to Grassmannian codebooks. 
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Fig. 1. The system model used in this letter. The source transmits in the first time slot, and the relay transmits in the second 
time slot. The relay is shown with separate transmit and receive antennas for convenience; this assumption is not made in the 
analysis. 
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Fig. 2. BER performance for several MIMO amplify-and-forward beamforming strategies. 
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Fig. 3. BER performance for uncoded BPSK versus Es/Nq for i.i.d. Rayleigh fading. At liigh SNR, thie slope of the curve 
approaches ~Ns{Nr + Nd) = —Nd{Ns + Nr) = —8, which, as shown in Section IV, is the full diversity order of the 
channel. 
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